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Abstract 

Approximate scattering functions for polydisperse ionic colloidal fluids are 
obtained by a corresponding states approach. This assumes that all pair cor- 
relation functions ^^/^(r) of a polydisperse fluid are conformal to those of an 
appropriate monodisperse binary fluid (reference system) and can be gener- 
ated from them by scaling transformations. The correspondence law extends 
to ionic fluids a scaling approximation (SA) successfully proposed for nonionic 
colloids in a recent paper. For the primitive model of charged hard spheres in 
a continuum solvent, the partial structure factors of the monodisperse binary 
reference system are evaluated by solving the Orstein-Zernike (OZ) integral 
equations coupled with an approximate closure. The SA is first tested within 
the mean spherical approximation (MSA) closure, which allows analytical so- 
lutions. The results are found in good overall agreement with exact MSA 
predictions up to relevant polidispersity. The SA is shown to be an improve- 
ment over the "decoupling approximation" extended to the ionic case. The 
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simplicity of the SA scheme allows its application also when the OZ equations 
can be solved only numerically. An example is then given by using the hyper- 
netted chain (HNC) closure. Shortcomings of the SA approach, its possible 
use in the analysis of experimental scattering data and other related points 

are also briefly addressed. 
05.20.Jj;61.10.Ex;61.12.Eq 
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I. INTRODUCTION 



Colloidal suspensions of charged particles represent a special class of ionic fluids |T]-Q. 
Unlike solutions of simple electrolytes such as NaCl, charged colloidal suspensions are highly 
asymmetric mixtures, containing both macroions and microions with large size and charge 
differences. Moreover, macroions often exhibit "polydispersity" , which means that particles 
of a same chemical species are not necessarily identical, because their size, charge or other 
properties may be spread over a large spectrum of values (chemical species whose particles 
are all identical are then referred as "monodisperse"). The presence of only one polydisperse 
macroion species is sufficient to make the colloidal suspension a mixture with a very large 
number p of components. The peculiar features of this "colloidal regime" , namely asymmetry 
and polydispersity, give rise to a variety of phenomena concerning microscopic ordering, 
phase behaviour, diffusion, and so on. 

Experimental information on the structure of such fluids can be obtained from Small 
Angle Scattering (SAS) techniques, by using light, neutrons or X-rays. However, when a 
signiflcant degree of polydispersity is present in the sample, the interpretation of experi- 
mental data for scattering intensity is hardly a simple task. In fact, polydispersity and 
large size-charge differences represent a serious challenge to the available theoretical tools. 
Monte Carlo or molecular dynamics simulations for polydisperse colloidal fluids involve very 
large numbers of particles. Moreover, large size asymmetries at high densities may cause 
ergodicity problems. On the other hand, integral equations (IE) of the liquid state theory 
are analytically solvable only under special conditions, whereas their numerical solution for 
mixtures with large numbers p of components, such as the polydisperse ones, would require 
large systems of non-linear equations. As a consequence, apart from very few peculiar cases 
1^, IE numerical studies on multicomponent fluids are usually restricted to p <^ 10. Finally, 
under the highly demanding conditions of colloidal suspensions non- convergence problems 
of the algorithms may often arise. 

The present paper will focus on the effects of polydispersity in SAS from ionic colloidal 
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mixtures, in the framework of IE theories based upon the Ornstein-Zernike (OZ) equations 
with approximate closures. Our study refers to the simplest polydisperse case with only two 
ionic species: monodisperse microions and macroions with both size and charge polydisper- 
sity. We shall refer to this system as polydisperse binary (two-species) ionic mixture. 

To overcome the impossibility of investigating polydisperse systems when lEs have to be 
solved numerically, one has to reduce the number of components and replace the study of a 
polydisperse fluid with that of a nearly equivalent but much simpler system. One possibility, 
not adopted in this paper, is to neglect microions altogether and approximate the fluid as 
a system of macroions interacting through a repulsive screened Coulomb potential, which 
implicitly takes into account the contribution of the neglected particles . A further refine- 
ment of this viewpoint is to build up an equivalent effective mixture with p' <^ p new 
components, whose molar fractions and diameters are determined by replacing a continuous 
distribution of macroion sizes with an appropriate p'-component histogram. Usually, p' = 3 
is already sufficient and therefore the problem is reduced to get a numerical solution of lEs 
for a three-component macroion mixture. This procedure could be easily extended to in- 
clude monodisperse microions and its counterpart would involve a four-component mixture 
(three for the macroions plus one for the microions). This method can be expected to be 
quite accurate, but it would demand a sizeable amount of numerical work. 

In this paper we present an even simpler approach, which requires the solution of only 
two-component lEs. We shall show that the problem of a polydisperse binary ionic mixture 
can be reduced to the study of a monodisperse binary ionic mixture, with microions and 
all identical macroions. The solution for such a reference system is the "starting" point for 
several approximations of increasing accuracy. Our main purpose is to show that, at the 
end of this hierarchy, accurate approximate scattering functions for a polydisperse binary 
mixture can be easily calculated with moderate numerical work upon using a corresponding 
states theory. Our method hinges on a conformality argument, which assumes that all 
pair correlation functions of the polydisperse fluid have essentially the same "shape" of their 
monodisperse binary counterparts and can be generated from them by means of simple 
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scaling transformations. This correspondence law is the extension to ionic mixtures of a 
scaling approximation (SA) successfully proposed for nonionic colloids in a recent paper 
1^. This is a non-trivial extension, since the good performance of the SA for the short- 
range potentials of nonionic colloidal fluids examined in Ref. ^ (uncharged hard sphere and 
Lennard- Jones interactions) does not automatically ensures the same success in the presence 
of long-range Coulomb attractions and repulsions. 

To properly treat both macroions and microions on the same footing, the colloidal sus- 
pension will be described by the primitive model (PM) of electrolyte solutions, which depicts 
all ions as charged hard spheres embedded in a dielectric continuum representing the sol- 
vent. The new SA will be tested, for the PM, against results from an analytic treatment of 
polydispersity, which is exact within the mean spherical approximation (MSA) closure for 
the OZ integral equations. In the PM-MSA case, the OZ equations were solved analytically 
many years ago |]8|-[TI|, and, more recently, a closed analytical formula was obtained for the 



scattering intensity from charged hard sphere fluids with any arbitrary number of compo- 
nents |[T^. An essential feature of the SA is that, because of its simplicity, this scheme can 
be applied equally well to combinations of potential models and closures for which only a 
numerical solution of lEs is possible. It is therefore possible for instance, as we shall ex- 
plicitly show, to couple the SA with the hypernetted chain (HNC) closure, which is more 
accurate than the MSA one for ionic fluids. In these cases the SA becomes a valuable new 
tool to predict properties of polydisperse colloidal suspensions in a very simple way. 

The paper is organized as follows. In the next Section the basic formalism of the small an- 
gle scattering and integral equation theory is briefly recalled along with the primitive model 
for polydisperse ionic fluid in Section III. In Section IV our corresponding states treatment 
of scattering functions is presented in detail, together with two simpler approximations. The 
exact MSA analytical expression for the scattering intensity from charged hard spheres is 
also reviewed and some of its predictions for polydisperse fluids will be reported. In Section 
V numerical results from the proposed approximations are compared in detail within the 
MSA. The performance of the SA with the HNC closure will then be addressed and few 



remarks will be included in the conclusive Section VI. 



II. SMALL ANGLE SCATTERING AND INTEGRAL EQUATION THEORY 
A. Scattering intensity and structure factors 

An ionic colloidal solution is formed by macroions and microions suspended in a homoge- 
neous solvent. Usually, this suspending fluid is formed by very small particles (with respect 
to the macroions) and is then modelled as a continuum, characterized by a given dielectric 
constant and an uniform density of scattering matter. 

According to the scattering theory, the intensity of the scattered radiation (light, neu- 
trons or X-rays) is proportional to the ensemble or time average, of |n(q)|^ over all possible 
equilibrium configurations of the sample particles. Here q is the exchanged wave vector and 
n{q) is the three-dimensional Fourier transform of n{r), a quantity related to the density of 
scattering matter at the position r inside the sample. For neutrons n(r) is the scattering 
length density: n(r) = bkS{r — r^), where bk is the scattering length of the k—th nucleus 
located at and 5 denotes the Dirac delta function. For X-rays n(r) coincides with the 
electron density, whereas for light it becomes the refractive index. In the continuum solvent, 
n{r) has an uniform value no. 

In addition to the continuum solvent hypothesis, we assume that inside each ion 
(macroion or microion, indifferently) the scattering matter has a well-defined boundary, i.e., 
there is a scattering core with a well-defined scattering volume, not necessarily coincident 
with the particle volume. While the former depends on the particle-radiation interaction, 
the latter is determined by the interparticle repulsions and is well-defined only in the pres- 
ence of hard body repulsions. The definition of a volume for particles with soft repulsions 
(e.g. Lennard- Jones particles) requires in fact some arbitrary and non-universal convention. 

In the case of suspended particles with spherically symmetric interactions (homogeneous 
and isotropic fiuid) and spherical homogeneous scattering cores, the SAS theory yields the 
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following expression for the total scattering intensity I{q) of a j9-component solution in a 
volume V 

p p 

R{q)^ I{q)/V = pY,Y.V^pPMFMSap{q)^ (1) 

a=l (3=1 

as a function of the magnitude of the exchanged wave vector, q = (47r/A) sin {6/2) , with A 
being the wavelength of the incident radiation and 6 the scattering angle. The Rayleigh ratio 
R{q) is the total scattering intensity per unit volume (also called the differential scattering 
cross section and often denoted by ^(q))- In Eq- (0) P is the total number density, while 
and Fi,{q) are the molar fraction and the form factor of species u respectively. Fi,{q) is 
related to the distribution n{r) of scattering matter inside particles of species u and we can 
express it as 

f.(9) = V;--'K-n,)ML^g^, (2) 

fjscatt ];5gj^g f^i^Q diameter of the scattering core, V'J*'^"** = | (a^^"**)^ its volume, rii, the uniform 
scattering density of species u and the difference Ui, — no its "contrast", while ji{x) = 
(sinx — a;cosa;)/x^ is the first-order spherical Bessel function. Finally, the functions Sapiq) 
are the Ashcroft-Langreth partial structure factors [|I^] 

Sccp{,q) = Sap + Py/XaXp hapiq), (3) 

where 6af3 is the Kronecker delta and hapiq) is the three-dimensional Fourier transform of 
the total correlation function hapir) = (7a/3(r) — 1. Here, gapi'r) is the radial distribution 
function (RDF) between two particles of species a and (3 aX a. distance r. 

In addition to the scattering intensity, it is then convenient to define a "measurable" 
structure factor |Q as 

p p 

Sniq) = XI XI V^»Xi3Waiq)w(3{q) Sapiq), (4) 

a=l 13=1 

with weights 



the brackets meaning = ^a^aF^ig)- The relationship between R{q) and Sm{(i) is 

Riq)=p{F'iq))SMiq). (6) 

From the theoretical point of view, we will obtain the partial structure factors Saisiq) by 
solving lEs for the /^^/^(r). 

B. Integral equations 

The OZ integral equations of the liquid state theory for ]9-component mixtures with 
spherically symmetric interparticle potentials are |p!4|,p!5| 



hap (r) = Caf3 {r) + ^Pu dr' c^u {r') /i^./? (|r - r'|) , (7) 



where the Cap {r) are the direct correlation functions and pu = x^p is the number density 
of species u. These equations can be solved only in combination with a further relationship 
between ha^ (r) and Cq,^ (r) . The formally exact expression of this "closure" is 

Cap (r) = exp [-{kBTr^UaiB (r) + -fap (r) + Bap (r)] -1-7^/3 (r) , (8) 

where Uap {r) is the interparticle potential, ks is Boltzmann's constant, T the absolute 
temperature, 'ja/s i^) = hap (r) — Ca/3 (r) and the "bridge" functions Ba/s (r) are functional 
of haf3 (r) and higher order correlation functions. In practice however, the exact Bap (r) 
cannot be calculated, and several approximations proposed for these functions define a 
corresponding series of approximate closures . 



The possibility of solving analytically the OZ equations depends on both the potential 
model Ua/s (r) and the chosen closure. Once that the lEs have been analytically or numeri- 
cally solved, the partial structure factors Sa/3 (q) can be obtained from Eq. (|^). 

III. PRIMITIVE MODEL FOR POLYDISPERSE IONIC FLUIDS 



We are interested in studying polydispersity effects by properly considering both 
macroions and microions on an equal footing. The simplest possibility is the primitive 
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{keT) (r) 



model (PM), well known in the theory of electrolyte solutions. It consists of an electroneu- 
tral mixture of p different components, represented by charged hard spheres embedded in a 
continuum solvent of dielectric constant e. The species a, with diameter cTq,, has molar frac- 
tion Xa and electric charge ZaC (e is the proton charge and Za the valency). The interparticle 
potential u^p (r) is defined by 

f 

+ 00 for r < = |((Tq, + (Tp) 

(9) 

ZaZpLs/r for r > a^p 

where Lb = e'^/leksT) is the Bjerrum length. The electroneutrality condition requires that 

The PM can also be used for polydisperse colloidal suspensions. In a "discrete repre- 
sentation" of polydispersity, a polydisperse two-species fluid is described by a p-component 
mixture {p ^ 1), in which the monodisperse microions (chemical species 1) are the first 
component, with diameter ai, charge zi and molar fraction xi, while the remaining p — 1 
components correspond to different varieties of the single macroion species (chemical species 
2). It is often convenient to adopt a "continuous representation" of polydispersity, with 
p ^ 00 and a continuous spectrum of values for the macroion "disperse" properties (size, 
charge, etc.). In such a continuous-mixture formalism, we assume that macroions have a con- 
tinuous distribution of diameters a around an average one, denoted by {(7)2 ■ For simplicity, 
we make the further reasonable assumption that the charge polydispersity of macroions is 
fully correlated to the size polydispersity. This can be easily accomplished by choosing 
the charge (or valency) of each macroion to be proportional to its surface area, i.e., 

-V 

.W2 



^2(a)=ZH J-^) , (10) 



where ^(o-)^ is the valency of the macroions having diameter (cr)2 . Both size and charge dis- 
tributions of the macroions are therefore governed by a single independent variable, namely 
the macroion diameter. The polydispersity of the macroions can then be expressed by a mo- 
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lar fraction density function, J92 (o") = x^2^f2{o'), where f2{o') is an appropriate distribution 
normalized to unit, while X2°* = 1 — xi is the "amplitude" of p2 (cr) . 

In the passage from a discrete to a continuous representation of polydispersity, the molar 
fractions Xu are replaced by dx = p2 (c) da, the fraction of macroions having diameter in the 
range (a, a + da), and the sums ^t^--- become integrals / da p2 (a) ... The average of a 
quantity Y over the macroion distribution is therefore written as 



-J^ 5^ ^ -L / da p2 (a) Y{a)= I da (a) Y{a), (11) 



where £2 denotes the set of indices corresponding to the macroion components. The average 
of Y over the whole set of suspended particles is then: (Y) = X]!^=i ^u^u = XiYi + x^2^ 0^)2 ■ 
For /2(o") = f{a; {a)2 , s), we use the Schulz or gamma distribution 

/(^;W,^) = r^^'^-'e-''^ («>1), (12) 



where F is the gamma function [|16|, while a = b = a/ (a) are related to the mean 



value {a) and the relative standard deviation s = y {a^) — (cr)^/ {a) , which measures the 
degree of polydispersity (0 < s < 1). The choice of the Schulz distribution is a popular one 
in colloidal theory because of its mathematical properties. For s —>■ 0, it reduces to a Dirac 
delta function centered at (a) (monodisperse limit). For small values of s, f{a) is similar 
to a Gaussian distribution, while for larger polydispersity it becomes rather skewed ^L7\ . 



Unlike the Gaussian function, the Schulz distribution is defined for positive values of a only. 
Moreover, this distribution allows a straightforward analytical evaluation of simple averages 
of the kind displayed in Eq. (|TTp. In particular, the first three moments of the distribution 
f2ia) are: {a)^ , {a'^)2 = (1 + s^) {a)l and {a^)^ = (1 + s^) (1 + 2^^) {a)l , while use of Eq. 
(p!0|) yields {z)2 = Z(a)^ (1 + s"^) ■ These analytical results can be conveniently inserted into 
the expressions for the electroneutrality and the packing fraction rj, i.e., 

x,zi + xf' {z)2 = 0, (13) 
r/ = (vr/6)p(xia? + xf (02)- (14) 
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The microion packing fraction is rji = (vr/6) picrf, while its macroion counterpart is 772 = 
(7r/6) P2°* (o"^)2 ' ^ith P2°* = p2;2°*- From Eq. (|13|) and 0:2°* = 1 — one then gets 



which shows that xi is fully determined by zi and {z)^ ( or equivalently zi, Z(^^),^ and s). 

A final remark is in order. In evaluating the averages of more complex quantities any 
analytical integration becomes a formidable or impossible task and numerical integration 
brings back to discrete expressions. For this reason, in the following we shall continue to 
employ the discrete notation under the implicit convention that Xa = a;2°'/2(o"a)Acr for the 
macroion molar fraction {Aa is the grid size in the numerical integration). 



Z(a)„ [l + S^ 



(15) 



IV. APPROXIMATIONS AND EXACT EXPRESSIONS 
A. Corresponding states and scaling approximation 

Interparticle potentials are said to be conformal when they have the same "shape" , and 
systems with conformal interactions are called conformal substances Analytically, the 
conformality of a set of potentials means that all their expressions can be generated from a 
single functional form by appropriate scaling of distances and potential parameters (particle 
sizes, energies, charges, etc.). 

The simplest example refers to pure fluids, when the potential Ua of any species a, in a set 
of substances, depends on only two parameters and can be written as Ua{r) = Sa u{r/aa) , 
where cTq, and Sa are a characteristic length and energy respectively, while m is a dimensionless 
function of the dimensionless ratio r/cTQ. The form of Uci{r) implies that all properties of that 
set of conformal fluids can be written in terms of dimensionless reduced variables, and it 
leads to the "corresponding states principle" commonly found in textbooks 0: all conformal 
pure fluids at the same dimensionless density and temperature have identical dimensionless 
pressure. The RDF of a pure fluid of species a in a group of conformal substances can be 
written as 
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ga{r;p,T;aa,ea) =g 




(16) 



where ^ is a universal function for such a group. If one among these fluids is arbitrarily chosen 
as reference system and its properties are labelled with the subscript 0, then its potential is 
uo{r) = eo M (r/cTo) and its RDF is given by go{r;p,T; ao,eo) =g{r/(jQ]pal, kBT/eo) . From 
Eq. (|16D one then gets 



where we have introduced dimensionless scaling factors = Ot^joa and = Ea/sQ. This 
result is tantamount to say that, if one knows the RDF of a reference fluid characterized 
by potential parameters cro,eo; then it is possible to derive the RDF of any conformal fluid 
of species a, with potential parameters cTq, The value of ga at r in a thermodynamic 
state (p, T) is equal to the value of go at the scaled distance A^r, in the corresponding state 
{p/Xa, T / with scaled density and temperature. For instance, if a^i > ctq and Ea > Eq, then 
the corresponding state has a greater density and a lower temperature. Using the definition 
of the potential of mean force, W = — ksT In g, Eq. ([T^ could also be cast in the form 



For pure fluids then conformality of the potentials implies conformality of the potentials of 
mean force and hence of the RDFs. The potential of mean force between two particles is the 
sum of the direct pair potential plus an indirect interaction, due to all the remaining fluid 
particles and averaged over all their possible equilibrium configurations. Finally, a similar 
property holds true for the structure factors as well 



The scaling correspondence in g-space is that the value of Sa at q is equal to the value of 5*0 
at A~^g (in a different thermodynamic state). 

On the other hand, for mixtures conformality of potentials does not necessarily ensure 
conformality of RDFs in the same simple way. Nevertheless corresponding states arguments 



ga{r]p,T;aa,Ea) = goiX^r ; p/Xl,T/^a]cro,£o), 



(17) 



Wa{r;p,T;aa,Ea) = Wo{Xar ; p/A^, T/^„; cxo, ^o)- 



(18) 



Saiq; p,T;aa,Ea) = So{Xjq ; p/ Xl,T /^a, ctq, Eq) . 



(19) 
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have often been exploited in the hquid state theory, to postulate approximate conformal- 
ity relations between mixture and pure RDFs p!^JT9|] . Only recently, however, this kind 
of approach has been applied to polydisperse fluids and a Scaling Approximation (SA) has 
been proposed for nonionic colloidal suspensions 0. In the SA theory is possible to ob- 
tain rather accurate structure factors for a "polydisperse one-species" fluid of uncharged 
spherical particles, by first evaluating the RDF qq of an appropriate "monodisperse one- 
species" (pure) reference fluid and then generating all the p{p + l)/2 different RDFs of 
the mixture by taking the values of the single at suitably scaled distances. The present 
work is aimed to extend this SA scheme to polydisperse ionic colloidal suspensions. It 
employs two-species fluids with both positive and negative ions, in order to satisfy the elec- 
troneutrality condition. As a reference system for the "polydisperse binary" fluid a suitable 
"monodisperse binary" (M2) mixture is required, where species 1 coincides with the mi- 
croions and has their density, size and charge, (pi™; -2^1") = (pi; -2^1); while the distri- 
bution of macroions is replaced by a single "average" component (species 2) with parameters 
(P2"; o"2™5 -2^2")- '^^^ choice of this reference fluid will be discussed later. Note that the set 
of parameters {p\''^, p'f'',T;af'',zf'',af'',zf'^) can be reduced to (p^'"^, T; ai, ^i, a^''^, 4''^), 
since x^*™ is automatically fixed by the charge ratio through the electroneutrality condition 
as: xf"" = (1 - zi/zf'^y^ . 

Our approximation consists in assuming that all RDFs of the polydisperse ionic mixture 
are conformal with the RDFs of the monodisperse binary fluid, which means that 

g^p (r; p, x,T; {a,s}, {z,s }) ^ g'Zm, {Kpr; p, T; a^,z^, af'^, z]^^) , (20) 

where x, {cr^^}, {2^5 } represent the complete set of molar fractions and potential parameters, 
pb- = p,Kp^ a'ZmJ^ap. With a'Z = « + <)/2, «,/? = !, and 

{1 when u = 1 
(21) 
2 when u E £2 

{£2 was already defined in Eq. ([TT])). The correspondence law ( PP] ) provides the recipe 
for generating all the p{p + l)/2 independent RDFs of the polydisperse fluid starting from 
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the three RDFs of the monodisperse binary mixture. It exphcitly reads gu (r) ^ (r) 
for microion-microion pairs, gi/3 (r) ^ g\™ (crj'™ r/ai/s^ , P E £2 for microions-macroions, 
gai3 (r) ~ (^2™ (cTj™ r/aais) , a,P E S2 for macroions-macroions. 

Our choice of Xap for scahng the distances imphes that, when r < cxc/j, one gets r'^p = 
XapT < cr^^^ and consequently ensures the correct hard core conditions, g^p ij) = for 
^ < o"q/3- The excluded volume effects, very important for the structure of condensed phases, 
are thus properly taken into account by the SA. 

Since the Fourier transform of h^^^^{Xaf3r) is h^^^{X~^q) , it is clear that S^l 



■ibin 



^mamfj + P°^^ \/ ^mZ^rri^ ^n^mp- Under the assumptiou that p*^™ = p and upon using Eqs. @ 



1^ ■-niamp 

and (|20|) one then obtains 



SA 



Sal3 + 



X: 



(22) 



where S^^^ (g) is a shorthand notation for S^^^^ (g; p, T; (Xi, Zi, 0-2™, z^™^ , which will be 
exploited hereafter unless otherwise specified. Eq. (|), within this approximation, takes the 
form 



+xT ^ ^ 



tot ^tot^c.iq)wpiq) ( ^ 



obin 
'-'22 



a/3 



(23) 



+ 2^/xiX^2^ Wi 



_.bin 
'-'12 



Eq. ( p3D is the basic result of the paper. It provides an expression for the measurable structure 
factor of the original polydisperse binary mixture, once that the partial structure factors of 
the reference monodisperse binary mixture are known. In the limit of vanishing charges and 
no microions it reduces to the one found in Ref . . The scattering intensity per unit volume 
R{q)^^ is then obtained by multiplying Spf{q)^^ by p {F^{q)) . 
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B. Choice of the monodisperse binary mixture 



As reference system, we select a monodisperse 2-component (M2) mixture which mimics 
the polydisperse p-component fluid. We assume that species 1 coincides with the microions 
and hence (p^'™, -Zi") = (pi, (Xi, zi), which implies the equality of the microion packing 
fraction, i.e., rj^^^ = rji. Then we replace the polydisperse macroion species, containing 
p — 1 components, with a monodisperse macroion species 2, containing a single "averaged" 
component. To determine its parameters (pg™, o"2™5 ^2^)5 "^^ require that 



^tot 

P2 5 



pizi + p^'°4'° = 0. 



2 ' 



(24) 



The first two equations guarantee that the total number of macroions and their packing 
fraction in the M2 mixture are the same as in the polydisperse fluid (^^2™ = V2)', the third 
one is the electroneutrality condition for M2. Combining Eq. ( p4|) with pY^ = pi, one finds 
the solution 





= P> 


< af- 


= (-)f 


^bin 
^2 


= (^)2 • 



and 



X 



bin 



(25) 



In this way, the definition of the set of M2 parameters (p^™, T; ai, Zi, cr2™, 22^^) is complete. 

Choices other than (^) are clearly possible. We have explicitly worked out few of them 
and found that they do not significantly alter the final numerical results. Eq. (EH) has then 



been privileged on the basis of its simplicity and natural physical interpretation. 

Besides being used as a reference system for SA, the M2 mixture may itself be regarded 
as the simplest approximation to the polydisperse p-component fluid. The corresponding 
measurable structure factor would then be 



SMiqr' = 1 + x.wliq) [S'^riq) - l] + xr[wY^{q)]' [ S^^ (q) - l] 

+2./^w,{q) wf-{q)S^^{q). 
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(26) 



which simply corresponds to approximate the original polydisperse binary mixture with a 
plain mono disperse binary mixture. 



C. Extended decoupling approximation 

To emphasize the role played by the scaling of distances in the SA, let us consider the 
simpler case of no scaling. This can be obtained from the SA expressions by setting X^p = 1 
everywhere. The result corresponds to an approximation which provides an exact evaluation 
of all form factors of the polydisperse system but assumes that the RDFs can be replaced 
by a set of only three effective RDFs of a monodisperse binary ionic fluid. Hence, in the 
previous language we have gu (r) ^ g\™ ( r) , gip (r) ^ g^^ (r) and ga/3 (r) ~ (^2™ (r) for 
micro-micro, micro-macro and macro-macro ionic pairs, respectively. Eq. (pSj) simplifies to 

SMigf"^ = 1 + x.wliq) [S^riQ) - 1] + i^ml [ (q) - 1] 

(27) 

+2./^w,{q) {w{q)), (g) . 

The superscript EDA means Extended Decoupling Approximation, since this approximation 
may be reckoned as an extension to polydisperse ionic colloids of the "decoupling approx- 



imation" (DA), proposed by Kotlarchyk and Chen pOf for nonionic fluids and well known 



to the small angle scattering experimentalists. The EDA may also be regarded as a special 



limiting case of the binary substitutional model proposed by Nagele et al. |^ for a different 
colloidal model with two polydisperse macroion species and no microions. 



D. MSA closure and analytic expressions 

The expressions we have previously derived for SA, EDA and M2 are clearly independent 
of the approximate "closure" chosen for solving the OZ equations. One then expects that 
an improvement in the selection of the closure would provide increasingly accurate results 
for the polydisperse colloidal suspension. In the present paper we shall focus mainly on 
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the MSA closure, to take advantage of its analytical properties. Another example will be 
considered in Section V-C. 

For the PM, the MSA consists in adding to the exact hard sphere condition, gap (r) = 
or haj3 (r) = — 1 when r < aaj3, the approximate relationship (closure) 

Cap{r) = -(ksTy'^Uapir) for r > aap, (28) 

which is asymptotically correct for r — oo. The advantage of the MSA closure is that 
the corresponding OZ equations for the PM were solved analytically some times ago [§. 
Senatore and Blum ||22| employed MSA expressions for the partial structure factors Sap{q) 



to calculate numerically Sm{(i) for charged hard spheres with either size polydispersity or 
charge polydispersity. More recently, a closedMSA formula for Sm{(i) by-passing the explicit 
calculation of the partial structure factors was obtained in Ref. |12[. This was the extension 



to ionic systems of an analogous expression for polydisperse uncharged hard spheres in the 
PY approximation ||23|. For the sake of completeness, the MSA analytic expression for the 
scattering intensity is reported in Appendix where some misprints appearing in Ref. 



are also corrected. The MSA closure yields analytic expressions for both 5"^°^^ (q) and 
Sm{(i) depending on a single screening parameter 2r, which in turn has to be determined 
self-consistently. 

A well known drawback of the MSA is that, for dilute solutions of highly charged parti- 
cles, it may predict unphysical negative values for g^p (r) near the contact distance aap or 
in a neighborhood of the first minimum. Some proposal have been advanced to heal this re- 



striction p4|-p6[| . For simplicity however, the emphasis of the present work will be mainly on 
concentrated suspensions of weakly charged particles. In this regime the MSA is reasonably 
accurate, with the Coulomb part of the potential being only a perturbation with respect to 
the hard sphere one. The above remark is nevertheless by no means a limit to our method 
which could be easily associated to more accurate closures such as the "hypernetted chain 
approximation" (HNC), corresponding to take Bap (r) = in Eq. (|^), or the self-consistent 
mixing scheme (HMSA) proposed by Zerah and Hansen for potentials with attractive terms 
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1^ , [T5[] . Clearly, in the HNC or HMSA integral equations the monodisperse binary reference 
fluid can be treated only numerically. This point will be further discussed in Section V-C 
where an example of such calculation will be provided for the HNC. 



V. NUMERICAL CALCULATIONS 

In order to display the behavior of the scattering functions under some typical polydis- 
perse conditions, we numerically reproduced a realistic experimental environment. 

The micro ions were given a valency zi = +1 and a diameter ai = 5 A (solvated counteri- 
ons), while we used an average macroion size ((7)2 = 100 A with relatively small charges ^(o-)^ 
in the range —20 ^0 (in e units) to ensure meaningful MSA results as previously discussed. 
We will increase this value up to Z(^^)^ = —50 later on using HNC. The scattering due to the 
microions is in principle not completely negligible and it might be also characterized by a 
different contrast with respect to the macroions. Nevertheless both contrasts were here fixed 
to the same value An = 4 x 10^° cm~^, which is typically found in neutron scattering from 
silica particles suspended in water In evaluating the form factors we further assumed 



af, = for all particles. A room temperature T = 298 K and the dielectric constant 
e = 78 of water result into a value Lb = 7.189 A for the Bjerrum length. All numerical 
calculations were performed for packing fractions rj = 0.1, 0.3 and polydispersity s = 0, 0.1, 
0.2 and 0.3 (the first value corresponding to the monodisperse binary mixture). We note 
that when s = 0, it is necessary to take |-22"^| < 30 for r] = 0.3 and |-22™| ^ 10 for i] = 0.1, 
to avoid unphysical negative values of the MSA ^fg" ^it contact. 

The three Schulz distributions, with polydispersity s = 0.1, 0.2 and 0.3, were discretized 
with a grid size Acr/ (a)^ = 0.02, and truncated at (Xcut/ (cr)2 = 1-56, 2.22 and 2.96, respec- 
tively. These cxcut values correspond to polydisperse mixtures with a number of macroion 
components equal to 79, 112 and 149, practically intractable with the available algorithms 
for solving lEs numerically. 
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A. Polydispersity and charge effects in exact MSA results 



Before analyzing the performance of the SA and other approximations, it is useful to recall 
how size and charges polydispersity affect the measurable scattering structure factor. This 
is achieved by using the closed analytical expression for Sm{(i) which is given in Appendix 
which is exact within the MSA. 

Fig. 1 depicts the effects of polydispersity on the measurable structure factor. Sm{(i)^^^ 
is plotted as a function of the dimensionless variable q {(j)^ for increasing values of s and 
fixed 7] = 0.3 and 2(a-)2 = —20. We note that as s increases at fixed t], p decreases. As 
expected, the effect of increasing polydispersity is three-fold: i) the oscillations on the tail 
of the curves are greatly reduced as a consequence of the destructive interference stemming 
from the several length scales involved; ii) the first peak is lowered, broadened and shifted 
to smaller q values corresponding to a larger typical distance between macroions-macroions 
nearest-pairs; iii) the q ^ limit is increased since highly dispersed particles can be more 
efficiently packed. All these effects parallel those observed in polydisperse nonionic fluids 
||7|JT7[] as well as in mixtures constituted of only macroions interacting through a repulsive 



screened Coulomb interactions ||2^, and they were already recorded even in the PM [ p^ . 

Next we check the effect of the charge. This is reported in Fig. 2, where the Spiiq)^^^ 
corresponding to r] = 0.3 , ^^o-)^ = —20, s = 0.3 is compared with that of the polydis- 
perse mixture of neutral hard spheres which results from "switching off" all charges and 
leaving all other parameters unchanged. As the charge increases, the main peak becomes 
higher and shifts to smaller q values, since its position is essentially determined by the 
macroion-macroion equilibrium distance which becomes larger in the presence of electro- 
static repulsions. The difference in the g — > behavior is also evident: the charges lower 
the S'm(?) values near the origin, as a consequence of the long-range nature of the Coulomb 
potential. 
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B. Scaling approximation plus MSA 



Our aim is now to display the performance of the SA when the partial structure factors 
^ruamfiil) of the reference M2 mixture are evaluated using the MSA closure. These SA-MSA 
results are compared with the exact MSA solution for polydisperse charged hard spheres 
previously discussed. We shall comment in the next subsection on a method we envisaged 
to avoid unnecessary repetitions in the IE calculations for the M2 mixture. 

In Fig. 3 the structure factor S'A/(g)^^~^^^ is shown for two different degrees of polydis- 
persity (s = 0.1 and 0.3) at low concentration and under weak charge conditions [r] = 0.1 
and = —10). The corresponding results from the M2-MSA and EDA-MSA approxima- 
tions are also reported for comparison. As expected, at small polydispersity (Fig. 3a) there 
are very little differences among all these curves, although the EDA yields a somewhat larger 
value for Sm{(1 = 0). This overestimation of the low-g scattering in the EDA becomes much 
larger as polydispersity increases (Fig. 3b). This is the same qualitative trend resulted in 
the DA for neutral systems |lT|,|3|J^,p9[] . At s = 0.3 it is apparent that the position of the 
first peak in the EDA follows that of the M2, whose maximum is shifted to larger q values 
with respect to the MSA result. On the other hand, the SA reproduces more accurately the 
position of the first peak and follows very closely the correct curve for q {a)^ ^ 6. 

The discrepancy in the low-g region, which is in fact the most interesting from the SAS 
point of view, can be more clearly seen in Fig. 4, where the scattering intensity per unit 
volume -R(g), calculated for the same parameters of Fig. 3 is displayed on a log- log scale. 
Nevertheless the SA performs overall rather well in all regions. 

In Fig. 5 the same functions of Fig. 3 are then plotted for a more concentrated suspension 
and higher macroion charges (r/ = 0.3 and zi^^j-^^ = —20), again for s = 0.1 and 0.3, while 
the corresponding results for R{q) are displayed in Fig. 6. Figs. 5a and 5b yield compelling 
evidence of the potentiality of the SA confirming the previous remarks. It is useful to analyse 
these results in the sequence M2 — > EDA SA. In Fig. 5a Sm{<i)^^^ exactly coincides with 
*S'm(?)^^ in the first peak region, but differs from it at the locations of the M2 minima and 



20 



at small q values. On the contrary, the SA curve is close to the exact MSA one everywhere. 
In Fig. 5b the M2 approximation largely disagrees with the MSA one. The EDA improves 
here on the M2 since it takes all form factors of the polydisperse fluid correctly into account. 
It exhibits a lower peak height and practically no subsequent oscillations. Nevertheless, 
the EDA has a dramatic low-g overestimate and it behaves poorly in essentially all regions. 
On the contrary, the SA is fairly accurate in the whole experimentally accessible g-range. 
Its performance at = 0.3 appears to be even more accurate then at r] = 0.1. All these 
features are quite remarkable if we recall that the SA, EDA and M2 curves have been 
obtained starting from the same partial structure factors S^^.^^{q). This fact clearly shows 
the crucial role played by the scaling of the distances. Its effect is to shift the first peak 
position to the right location and to dump all oscillations after the first peak. Physically 
it confirm the soundness of our "conformality" hypothesis as expressed by Eq. ( pOD and it 
shows that overlooking differences among macroion-macroion RDFs at contact (as it is done 
in the SA) is a reasonable assumption. 



C. Scaling approximation plus HNC 

To illustrate the possibility of applying the SA scheme even when the OZ equations 
admit only numerical solutions, we investigated SA with the HNC closure and analyzed two 
cases: r] = 0.1, s = 0.3 with zi^o-)^ = —10 and Z(^„)^ = —50, respectively (all other parameters 
were fixed as before). While the first case was already studied with the SA-MSA, the second 
one represents a situation, of low concentration and high charges, in which the MSA yields 
negative values of the macroion-macroion RDF at contact and hence cannot be utilized. 
This drawback is avoided by using the HNC closure. 

In both cases we solved the HNC equations for the parameter values of the corresponding 
M2 reference mixture, using an r-space grid size Ar/ {a)^ = 0.02 and a number of grid points 
J\f = 4096. This choice implies that q^^^ {(t)^ = 50n, with a small enough grid size in g-space, 
Aq = gmax/A/", allowing the implementation of a "trick" proposed in Ref. . In fact, the 
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expression for Sm{(i)^^, Eq. (p3D, would require, at each q, the evaluation of one term 
5j'}°(g), p-1 terms S^^"" {q(yip/o-^t) and (p- l)p/2 terms ^2^^° (ga^/j/fx^'"") (recall that p - 1 
is the number of macroion components). These cumbersome repeated calculations can be 
avoided. We calculated S^f", S"^^"" and S^^ at the grid points qi = iAq {i = 0, . . . ,Af - 1) 
only once, storing all these values in arrays. Although the grid points do not exhaust 
the whole set of qaif^/a^^ and gcrQ,/3/(T2™ values required in Eq. (^), the stored structure 
factors represent a fine sampling of these continuos functions. Therefore, if the Aq is small 
enough, the value of S"}*" (or 6*22°) at a certain point can be approximated with that at the 
nearest grid point with a negligible error. In this way the sums of Eq. (^) can be quickly 
performed. 

Figure 7 shows the Sm{(i)^^~™^ curves, along with the SA-MSA one for Z(^„)^ = —10. 
As expected, in the lower charge case the SA-HNC prediction is very close to the SA-MSA 
one, with only a slight shift in the first peak position. The SA-HNC structure factor with 
Z{a)^ = ^50 is a not trivial result: it refers to a polydisperse colloidal system with mean size 
ratio 0"! : {cr)^ =1:20 and mean (absolute) charge ratio \zi\ : |2(<j)2| =1:50. Unfortunately, 
in this case we cannot make a comparison with "exact" data. On the other hand, the lack 
of these data and the difficulty of generating them in a very asymmetric regime is just the 
strongest motivation for introducing approximate theories such as the SA. 

VI. CONCLUDING REMARKS 

In this paper the problem of computing scattering functions for polydisperse ionic col- 
loidal fluids has been addressed by integral equation methods. In the framework of the 
primitive model we have shown that, despite the complexity of these systems, surprisingly 
accurate predictions can be obtained with a limited numerical effort. We have successfully 
extended the scaling approximation introduced in Ref. for polydisperse fluids of neutral 
particles. The SA still works well when Coulombic (both repulsive and attractive) interac- 
tions are present, notwithstanding the strong charge-size asymmetries of the polydisperse 
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colloidal regime. Only the study of an appropriate monodisperse binary mixture (the M2 
reference system) is required for a complete characterization of the polydisperse system. 

Our corresponding states-like theory is based on the simple physical idea of conformality 
of all RDFs in the polydisperse mixture. All partial structure factors are generated by 
scaling their three counterparts of the M2 fluid. In the liquid state theory similar ideas have 
been widely exploited in the past [l^JT9(l but, to our knowledge, Ref. and this paper are 
the first application to polydisperse fluids. 

Clearly, the SA theory is accurate only in the average. In fact, the scaling is hardly 
accurate for each individual pair correlation ^'^/^(r). In particular, it incorrectly assumes 
the equality of the RDFs at contact for all macroion-macroion pairs as well as for the 
microion-macroion ones. However, an essential feature of our SA is that it correctly ensures 
Qapij") = inside the hard cores. These excluded volume conditions are crucial, as it is 
shown by the failure of the "extended decoupling approximation" which neglects them. The 
structure factors Sj^f{q)^^ turn out to be accurate in the first peak region and beyond; some 
inaccuracy, due to the harsh approximations of our theory, is found at low q values. Since 
the q ^ limit is related to thermodynamics, this means that the SA can be meaningfully 
exploited to extract structural but not thermodynamical predictions. 

Because of its simplicity, the SA can be safely employed also when the OZ integral equa- 
tions have to be solved numerically since its application to both different closures and dif- 
ferent potential for ionic colloids is feasible. These features are indicative that SA is a useful 
theoretical tool to investigate, to first approximation, the structure of polydisperse (nonionic 
and ionic) colloids under highly demanding conditions. The existence of a good approxima- 
tion which reduces the study of polydisperse fluids to that of an effective monodisperse one 
should not be underestimated. Real life colloids are always polydisperse to a certain degree 
and polydispersity always represents a challenge in the interpretation of experimental data. 
We hope that the SA will result particularly useful in the analysis of Small Angle Scatter- 
ing data, since it considerably outperforms the "decoupling approximation" , popular in this 
context, at the cost of a minimal additional effort. 
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It would be interesting to compare our theory with the approach proposed by D'Aguanno 
and Klein ^ . As already mentioned in the Introduction, these authors followed a different 
point of view and replaced the continuous Schulz distribution with an histogram containing 
a finite number of well-chosen diameters, thus reducing the polydisperse fluid to an effective 
mixture with a very small number p' of components. In this case however we expect a 
non-trivial increase in the numerical effort involved. In fact, to avoid the rapid increase in 
computational cost with increasing p' in the D'Aguanno-Klein mixtures, that approach was 
recently modified by Lado and coworkers by adding an orthogonal polynomial expansion 
technique , and afterwards by merging this with a thermodynamic perturbation scheme 



11]. A comparison of the SA with these alternative theories is left to future work. 



ACKNOWLEDGMENTS 

Financial support from the Italian MURST (Ministero dell'Universita e della Ricerca 
Scientifica e Tecnologica) through the INFM (Istituto Nazionale di Fisica della Materia) 
is gratefully acknowledged. One of us (DG) thanks Giorgio Pastore (Trieste) and Enrique 
Lomba (Madrid) for providing their HNC computer codes. 

APPENDIX A: MSA EXPRESSION FOR THE SCATTERING INTENSITY 

In this appendix we report the basic formulas involved in the MSA calculation of the 
scattering intensity from charged hard sphere fluids, as described in Let us introduce 
the following short- hand notations: 

{y}o^f;p,n = p(F), (Al) 

i/=l 

{F} = X^p.ne^^''=p(ye^^>, (A2) 

where X,y = qay/2. The MSA analytical solution depends on the screening parameter 2r, 
which must be determined numerically by solving the consistency equation 
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where 



m' = ^^Ls±pJ '\f;.f^ )\ (A3) 

L>=1 ^ ^ / 



Q I 1 + To- 



vr I 







n = A + -<; :;— }> , (A5) 



2 U + Tcr 

with A = 1 — r/. These quantities are also required to compute 



(A6) 



r l + Fa, 

In the hmit of point ions (all 0), 2r becomes the Debye inverse shielding length Kd 

of the Debye-Hiickel theory for electrolyte solutions, while for finite size ions it is always a 
lower bound (i.e. 2r < kd)- We also need 

A=^. (A8) 

= A^a.i,, , (A9) 

q J-iB 

^i'^ = A.e-^^^ (AlO) 

7. = 7f + 7^^\ (All) 

where = 3q{Xi) and = ?>ji{Xi,)/ X,^, with ^'0(2^) = sinx/x and ji{x) — {cosx — 
xcosx)/x^ being Bessel functions. 

The final expression for the scattering intensity per unit of volume is 

R{q)^Ri{q) + R2{q), (A12) 
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where 



+ 2Re [ Ci + {Fpj^ ca + {ap}^ Cicl 



(A13) 



R2 = {Iqifjo Icsl^ + 2Re [ {g7^}o cs + {qiajo c^cl + {g7/5}o C3C2 



(AM) 



Here, is the form factor given by Eq. Re [ . . .] the real part of a complex number and 
the asterisk denotes complex conjugation. Other necessary quantities appearing in these 
equations are 



to 



Cl 



C2 



C3 



(A15) 



ti ti ti 

(m = 1, . . . , 4) being the cof actor of the (1, m)th element of the first row in the following 
determinant 





1/2 
Pu a„ 


1/2/3 
pu Pu 


1/2 
pu qiu 


{F} 


1+ {a} 


m - 36/A + zg/2 




{aF} 


{era} 


1 + {aP} 




{AF} 


{Aa} 


{A(3} 





(A16) 



and where ^2 = (tt/O) {ct^Iq • 

Eqs.( [A4|) and ( [A16| ) correct the misprints appearing in the corresponding equations of 
Ref. [0. In the expression of given by Eq. (48) of that paper the factor tt was omitted. 
We also note that our definition of and f2 differs from Blum's original one [R] . 
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FIGURES 

FIG. 1. Polydispersity effects. Exact MSA structure factor Sm{q) of polydisperse charged 
hard spheres, at fixed packing fraction t] = 0.3, for different degrees of polydispersity s, (s = 
corresponds to the monodisperse binary case). Other parameters: cri = 5 A, zi = +1 (e units) for 
microions; ((7)2 = 100 A, zi^^^^^ = —20 for macroions. 

FIG. 2. Charge effects. The exact MSA structure factor Sm{(1) for polydisperse charged hard 
spheres with rj = 0.3, Z(^^^^ = —20, s = 0.3 (other parameters as in Figure 1) is compared with its 
exact PY counterpart for the corresponding polydisperse mixture of neutral hard spheres. 

FIG. 3. (a) Structure factor SM{q) for rj = 0.1, Z(^^)^ = — 10, s = 0.1 (other parameters as in 
Figure 1). Comparison of M2-MSA, EDA-MSA, SA-MSA and exact MSA results, (b) Same as 
(a), but for s = 0.3. 

FIG. 4. Scattering intensity per unit volume, R{q), using a log-log scale. Comparison of 
M2-MSA, EDA-MSA, SA-MSA and exact MSA results. The systems are the same as in Figure 3: 
(a) ri = 0.1, Z(^^)^ = —10, s = 0.1. (b) Same as (a), but for s = 0.3. 

FIG. 5. (a) Structure factor Sm{q) for rj = 0.3, z^^^^^ = —20, s = 0.1 (other parameters as in 
Figure 1). Comparison of M2-MSA, EDA-MSA, SA-MSA and exact MSA results, (b) Same as 
(a), but for s = 0.3. 

FIG. 6. Scattering intensity per unit volume, R{q), using a log-log scale. Comparison of 
M2-MSA, EDA-MSA, SA-MSA and exact MSA results. The systems are the same as in Figure 5: 
(a) rj = 0.3, Z(^^^^ = —20, s = 0.1. (b) Same as (a), but for s = 0.3. 

FIG. 7. SA-HNC predictions for the structure factor Sm{q) at rj = 0.1, s = 0.3, with 
^{'^)2 ~ ^^'^ = (other parameters as in Figure 1). In the zi^u^^^ = —10 case the 

corresponding SA-MSA curve is also plotted for comparison. 
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